quantitative analysis of the critical current due to vortex pinning by surface 

corrugation. 
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The transport critical current of a Niobium (Nb) thick film has been measured for a large range of 
magnetic field. Its value and variation are quantitatively described in the framework of the pinning 
of vortices due to boundary conditions at the rough surface, with a contact angle well explained by 
the spectral analysis of the surface roughness. Increasing the surface roughness using a Focused Ion 
Beam results also in an increase of the superficial critical current. 

PACS numbers: 74.60.Ec, 74.60.Ge, 61. 12. Ex, 74.70.Ad. 



INTRODUCTION 

The understanding of the nature of the critical cur- 
rent in superconductors has triggered a lot of work for 
many years 0,0- Practically speaking, quantitative pre- 
dictions of the critical current values is a key point for 
most of the applications. More fundamentally, the ques- 
tion of the interaction of the vortex lattice with pinning 
centers is also a prototype system for understanding the 
behaviour of an elastic medium submitted to a disorder 
potential 3] . It is well known that a perfect vortex lattice 
( VL) submitted to a bulk transport current is subject to a 
flow, that leads to a finite electro chemical field and gives 
rise to dissipation . The ability of the VL not to move 
when it is submitted to a transport current is generally 
explained by the pinning interactions between pinning 
centers and VL. Unfortunately, the exact nature of the 
pertinent pinning centers, and the way they are acting, is 
not straightforward. In soft samples where superconduct- 
ing parameters vary slowly in the sample, cristallographic 
bulk inhomogeneities are usually supposed to play the 
role of pinning centers Q{. One can notice that this col- 
lective bulk pinning description leads to quite complex 
critical current expressions and to a corresponding lack of 
quantitative interpretation and prediction. On the other 
hand, experimental observation of the strong influence 
of surface quality on VL pinning has been evidenced in 
pioneering experiments |5|. This surface influence can 
be described in the framework of surface roughness in- 
teracting with VL, as proposed by Mathieu-Simon (MS) 
in a continuum approach. For any real sample, which 
presents a surface roughness at the VL scale, the respect 
of the boundary conditions of terminating vortices at 
the surface imposes local bending. This defines a contact 
angle (8—0 for a flat surface) and leads to a near-surface 
supercurrent. At the sample scale, this offers a lot of 
metastable equilibrium states. The ability of the system 



to sustain a macroscopic supercurrent is then directly 
linked to the surface roughness via an average contact 
angle C . The macroscopic critical current (per unit of 
width) is simply given by 



i c (A/m) = e.sin(6 c ) 



(1) 



where e is the vortex potential, i.e. the thermody- 
namical potential describing the equilibrium state. "Pin- 
ning" is then nothing else than the consequence of vor- 
tex boundary conditions applied to a real surface 
This should a priori act in every sample. The remain- 
ing question is: what is the contribution of this surface 
critical current with respect to the overall critical cur- 
rent? The vortex potential can be measured from the 
reversible magnetic moment M = -Le <5V of a super- 
conducting volume V, or calculated using, for example, 
Abrikosov's solution close to B C 2 . Using an iteration pro- 
cedure proposed by E.H. Brandt |ic||. e can be also ob- 
tained over the whole mixed state via numerical solutions 
of the Ginzburg-Landau equations. Now if one assumes 
a physical value of C = 0-20 deg, a good order of mag- 
nitude of i c is deduced |3 ■ A careful inspection of the 
surface state should in principle enable to extract a more 
precise value of this critical angle. If the usual measure- 
ment of surface roughness is the rms (root-mean-square) 
height h of the surface bumps, the pertinent parameter 
is here more the distribution 9(x) = arctan (^) of local 
slopes over the width of the surface. The rms value of 
this angle is given by the integration of the spectral den- 
sity 2 = J k max See dk over the appropriate k boundaries 
See is the Fourier transform of the autocorrelation 
of 0(x) (Wiener-Khintchine theorem), which represents 
the spatial distribution of 6. The aim of the present 

study is to measure the critical current of a Niobium 
thick film, with various surface corrugation, and to com- 
pare the experimental values to those obtained with the 
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simple expression (1). e will be calculated using Brandt's 
approach, and 9 C compared to y/fjP. 

EXPERIMENTAL 

The sample used is a film of Niobium (thickness^ 3000 
A) deposited at 780C on a sapphire substrate by the ion 
beam technique. The film has a resistivity of about 0.5 
/ifi.cm at the critical temperature T c — 9.15 K and ex- 
hibits a low surface rms roughness (R a 5 nm), mea- 
sured by Atomic Force Microscopy (Nanoscope III, Digi- 
tal Instruments). Microbridges of W=10 /mi x L=30 /im 
have been patterned using a scanning electronic micro- 
scope, this irradiation step being followed by a reactive 
ion etching process. The critical currents have been mea- 
sured by mean of the standard four-probe technique, at 
the following temperatures of 4.2, 5.2, 6.2 K in the whole 
range of field covering the mixed state. The critical cur- 
rent values I c were determined with a voltage criteria of 

10 nV. 

RESULTS AND DISCUSSION 

Let us first discuss the general behavior of the I c {uf) 
data (fig. 1) for the virgin microbridge at different tem- 
peratures (lu is the vortex field inside the sample and 
corresponds to the vortex density n=^-) . Note that the 
demagnetization factor due to the geometry of the thin 
film renormalizes the apparent first critical field B c \ up 
to about 0.015 £? c ii>«tt- It implies that the mixed state 
is created for the lowest field value applied (m 30Gauss). 
Except for this peculiarity , I c (uj) exhibits the same field 
variation as the reversible magnetization curve of a type 

11 superconductor, in agreement with the expected varia- 
tion of superficial currents i c (equation 1 with 9 C w cte) . 
We also notice that the curves taken at 4.2, 5.2 or 6.2 
K are self similar, that is they can be superimposed by 
a mere rescalin g. As first noted in pioneering work on 
vortex pinning [12j, the change of critical current with 
temperature can be totally attributed to the change in 
primary superconducting properties. It is evidenced in 
figure 2 where the low field value of I c is shown to be sim- 
ply proportional to B c2 for the three temperatures and 
both for the virgin and for the damaged sample. This 
shows that the variation of I c with temperature is sim- 
ply due to the variation of the vortex potential e(T)(i.e. 
the variation of superfluid density), without the need of 
involving other thermal effects such as vortex thermal dif- 
fusion. Now if we want to verify quantitatively equation 
(1), we first need to know the vortex potential e. It is usu- 
ally approximated using the Abrikosov calculations from 
B C 2 down to 0.4-0.5 B C 2, and by the London expression at 
very low fields B > B c \. For low kappa superconductors 
such as pure Niobium (/t~l for a sample with properties 



very close to ours [l3j[14|), only the Abrikosov expression 
is quantitatively correct and it is necessary to use numer- 
ical calculations to solve the Ginzburg Landau equations 
over the whole field range, following for example Brandt's 
iterative method [l5j . The result is presented in figure 3, 
for k — 1. We can now deduce the critical angle needed 
to account for the measured critical current. For a pure 
surface pinning and following equation (1), the expected 
value of 9 C is given by 9 C w arcsin (i c /e). Using the nu- 
merical values of e, we find 9 C w 0.4-1 deg (see figure 
4), in agreement with the mere expectation of a physi- 
cal angle. The order of magnitude is promising, but the 
complexity of the (multiscale) surface disorder needs a 
careful surface analysis. We have therefore measured the 
microbridge roughness using AFM in tapping mode (see 
ng. 5). 

Following the simple analysis described in the intro- 
duction, we obtain the spectral density Sgg with the use 
of the Fourier transform of the autocorrelation of 6(x) 
(Wiener-Khintchine theorem). The main value of the 
statistically representative angle is given by the integra- 
tion of the spectral density E2: 9 2 = $ max S eg dk. The 
boundary of the integral have been chosen as the natural 
scale for the vortex lattice, considering that a vortex does 
not see a roughness less than its core size (k max = 
Note that the choice of this cutoff frequency does not 
significantly change the results. This calculation gives 
\/~6~ 2 ~ w 0.60 ± 0.18 deg. The agreement with the value 
deduced from the critical current measurements (a mean 
value of 0.70 ± 0.15 deg) is quite promising. More pre- 
cisely, this is, within error bars, what we obtained for 
the magnetic field values higher than about 1000 Gauss. 
It is worth noting that the main value that we calculate 
is statistically representative but gives also a value that 
is supposed to be independent of the frequency. We are 
fully aware that a more rigorous analysis should take into 
account a kind of matching effect between the VL period- 
icity and the scale of surface disorder. It is even possible 
to expect a peculiar variation of i c (u>) in the case of a 
very rough surface at a restricted spatial scale. In this 
respect, one can see that the highest angles observed for 
low fields B < 0.1 T (ao > 0.5 /im) are quite consistent 
with the highest angles observed in the surface profile for 
a periodicity of about 1 /im. 

We decided also to compare this virgin microbridge 
with one whose surface structure was modified. The idea 
was to use a Focused Ion Beam to etch its surface fol- 
lowing a controlled geometry. The expected shape was 
that of "corrugated iron" with 12 /im by 0.1 /im trenches 
regularly spaced by 1 /im. Also, the etched depth should 
be high compared to the initial roughness of the surface, 
that is 30nm here. Several attempts have been performed 
using different Ga ion doses, which should be as low as 
possible in order to minimize the effect of Ga irradiation, 
leading to the best control of the etched surface. Thus, 
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the final procedure was eight identical patterns (12/im 
x 0.10/xm x 0.03/mi), etched parallel using an ion cur- 
rent of 4 pA corresponding to a total ion dose about 150 
pC//j,m 2 . The sample was tilted by 45 degrees and the 
magnification used was 25kX. Figure 6 is a SEM image of 
the resulting etched microbridge; the trenches are evenly 
spaced and they yield geometric parameters close to those 
expected (width about 0.15 /im). Such low energy irra- 
diation leads also to an implantation of Ga + ions, but 
simulations using Monte-Carlo calculations [ljS] indicate 
that it affects only a range of no more than 100 A. Fur- 
thermore, we observe neither any change in the critical 
temperature (T c = 9.15-ftT), nor in the normal state re- 
sistivity (p„ w 0.bO[iQ.cm) and critical fields within ex- 
perimental accuracy. In order to evaluate more precisely 
the influence of the Ga irradiation, another microbridge 
was etched using a single rectangular pattern (12m x 8m 
x 0.03m), that is covering the whole width of the mi- 
crobridge. The sample was tilted of 45 and the total ion 
dose close to 190pC//mi 2 . Critical current and resistivity 
measurements were performed on this microbridge and 
were compared to those on the virgin microbridge. We 
observed the same properties and specially that the crit- 
ical current is very similar (within few percents) in the 
two cases. This confirmed that the Ga irradiation had 
no influence on the bulk physical properties of the film. 
The FIB treatment, contrary to highest energy irradia- 
tion, has not modified the bulk crystalline lattice of the 
material and that the important modification is in the 
surface structure. We can therefore use the term "sur- 
face" damages. The main obvious result is the increase 
of the critical current, as shown in fig. 7. Following the 
same procedure that we performed for the virgin micro- 
bridge, one obtains w 2.2 ± 0.3 deg in good enough 
agreement with the 1.5 ± 0.2 deg deduced from the crit- 
ical current values. We note that the kind of treatment 
we performed leads to an increase of the roughness for a 
periodicity of about 0.1-1 /im. In the inset of figure 7, one 
can evidence that the critical angle is almost unchanged 
for the highest magnetic field values u> > 3000 G. For 
those values, the intervortex distance is less than 0.1 /im 
and we find using the AFM that the treated microbridge 
exhibits the same kind of roughness as the virgin micro- 
bridge for this periodicity. Again, the spatial dependence 
of the surface roughness is certainly linked with the ex- 
act 9 c (u) variation, with u> fixing the spatial scale of the 
pertinent surface disorder. More work is needed to fully 
describe this problem. 

Finally, we conclude from this analysis that the critical 
current of our sample is given by equation (1) with a good 
agreement, for the virgin and for the surface damaged 
microbridges. We conclude that taking into account the 
surface defects as main sources of pinning enables to ex- 
plain the experimental critical current values. Expressed 
in the form of a surface critical current and due to the 



small surface corrugation of the Nb Film (a roughness of 
few nm rms), the critical current appears to have stan- 
dard and even relatively small values for this low kappa 
superconductor (i c w 10-30 A/cm at low fields). If one 
expresses the critical current in the form of a density 
as it is usually made, this leads to a high value (J c « 
0.5 — 1.5 10 6 j4/cm 2 ). We emphasize that this notion of 
density is by definition not justified in the case of a cur- 
rent flowing under the surface and not uniformly in the 
bulk. As an example, Niobium crystals with a thickness 
of 7.6 10 -3 cm exhibit a critical current density of roughly 
3.10 3 A/cm 2 at 0.1 T and 4.2 K 17]. This gives i c =13.1 
A/cm compared to 12.5 A/cm (our thin film) under the 
same (w,T) conditions. So, The surface critical current is 
almost the same. It follows that the difference of thick- 
ness makes this apparent (but not physically significant 
in terms of pinning "force") difference in the critical cur- 
rent density values. Note that the same remark apply to 
other types of clean superconductors [l8j |. 

Now, if one increases a lot the number of bulk defects to 
obtain a spacing say less than the intervortex distance, we 
can obtain a so-called hard superconductor. In this case, 
a bulk subcritical current can flow by a percolating-like 
behavior. Anyway, in this case, the critical current den- 
sity was shown, in a lot of cases, to be proportional both 
to the difference in equilibrium magnetization across the 
internal boundaries and to the area of internal interface 
per unit volume 0, |2£j, El ■ It is worth noting that this 
is the same kind of pinning by surface interactions, but 
in this case with internal interfaces. 

In summary, we have observed that the value of the 
critical current of a thick film of a standard type II su- 
perconductor is quantitatively explained with the vortex 
pinning by surface roughness. We have also observed the 
increase of this critical current caused by an increase of 
the surface corrugation. Furthermore, it gives a simple 
explanation for the high critical current density observed 
in this kind of clean thin films, compared to the mod- 
erated one observed in (thick) bulk crystals. We hope 
that it gives also evidence that the interaction between 
the surface corrugation and the vortex elasticity is a key 
point for the understanding of vortex lattice pinning and 
dynamics. 
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FIG. 1: The critical current of the virgin niobium microbridge 
as function of the vortex field for three different temperatures. 

FIG. 2: The critical current at 30 Gauss as function of the 
second critical field for the three temperatures and the virgin 
and surface damaged microbridges. 



FIG. 3: The variation of the vortex potential (or equilibrium 
magnetization) as function of the magnetic field for re=l cal- 
culated with Ginzburg-Landau equations and following the 
method explained in 10]. The two are respectively normal- 
ized over Hc2 and Bc2- The dashed line is the Abrikosov 
line. 



FIG. 4: The variation of the critical angle 6 C deduced from the 
arcsin(^). The straight line corresponds to the main value 
extracted from the surface analysis. 
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FIG. 5: AFM picture of the surface roughness for the virgin 
microbridge. In the inset is shown the corresponding variation 
of local slopes 8 over the width of the sample. 



FIG. 6: a/ SEM image of the treated microbridge showing 
the eight trenches (arrowed) in the niobium film, b/ Detail 
of one trench showing its geometric characteristics. 



FIG. 7: comparison of the critical currents for the virgin mi- 
crobridge and for the one degraded with the FIB. In the inset 
is shown the critical angle variation deduced from equation 
(1). 
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